m 



PROCEEDINGS 

°F SCIENCE 



m 
o 

(N 

a 



Polyakov loop analysis with Dirac-mode expansion 



> 

in 
in 
oo 
(N 

O 

cn 



Takumi Iritani*, Shinya Gongyo, and Hideo Suganuma 

Department of Physics, Kyoto University, Kitashirakawaoiwake, Sakyo, Kyoto 606-8502, Japan 
E-mail: 



iritani@ruby . scphys . kyoto-u . ac . jp 



In order to investigate the direct relation between confinement and chiral symmetry breaking in 
QCD, we investigate the Polyakov loop in terms of the Dirac eigenmodes in both confined and 
deconfined phases. Using the Dirac-mode expansion method in SU(3) lattice QCD, we analyze 
the contribution of low-lying and higher Dirac-modes to the Polyakov loop, respectively. In 
the confined phase below , after removing low-lying Dirac-modes, the chiral condensate {qq) 
is largely reduced, however, the Polyakov loop remains almost zero and Z3-center symmetry 
is unbroken. These results indicate that the system is still in the confined phase without low- 
lying Dirac-modes. By higher Dirac-modes cut, the Polyakov loop also remains almost zero 
below Tc- We also analyze the Polyakov loop in the deconfined phase above Tc- We find that the 
Polyakov loop and Z3 -symmetry behavior are insensitive to low-lying and higher Dirac-modes in 
both confined and deconfined phases. 
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1. Introduction 

Quantum chromodynamics (QCD) is the fundamental theory of the strong interaction, how- 
ever, its non-perturbative properties such as confinement and chiral symmetry breaking are not yet 
well understood. In particular, to clarify the correspondence between confinement and chiral sym- 
metry breaking is one of difficult and interesting subjects ^, ||, ^ ^ ^, ^. As an evidence of the 
close relation between them, lattice QCD calculation shows that simultaneous deconfinement and 
chiral phase transition at finite temperature [^. 

As shown in the Banks-Casher relation [^], the chiral condensate (qq) is proportional to the 
Dirac zero-mode density as 

(^^) = -lim lim7r(p(0)), (1.1) 

where p (A ) is the Dirac spectral density. Thus, the low-lying Dirac eigenmodes directly relate to 
chiral symmetry breaking, however, their relation to confinement is still unclear. 

Therefore, it is interesting to analyze confinement in terms of the relevant degrees of freedom 
for chiral symmetry breaking, i.e., the low-lying Dirac eigenmodes. For example, based on the 
Gattringer's formula [|l|, the Polyakov loop can be investigated by the sum of Dirac spectra with 
twisted boundary condition on lattice |^ ^, §]. In our previous studies [js), §], we formulated the 
Dirac-mode expansion method in lattice QCD, and investigated the Dirac-mode dependence of the 
Wilson loop and the interquark potential. It is also reported that hadrons still exist as bound states 
without "chiral symmetry breaking" by removing low-lying Dirac-modes |[T^, H |. 



In this paper, we investigate the Dirac-mode dependence of the Polyakov loop in both confined 
and deconfined phases at finite temperature, using Dirac-mode expansion method in lattice QCD 
[§, |6|, ^. In Sec. 2, we review the Dirac-mode expansion method, and formulate the Dirac-mode 
projected Polyakov loop. In Sec. 3, we perform the lattice QCD calculations for the Polyakov loop 
with Dirac-mode projection. Section 4 is devoted for the summary. 

2. Dirac-mode expansion method in lattice QCD 

Here, we introduce the Dirac-mode expansion technique in lattice QCD [||, ^, ^, and formu- 
lation of the Dirac-mode projected Polyakov loop. 

2.1 Dirac-mode expansion in lattice QCD 

Using the link- variable € SU(A'^r), the Dirac operator D = y^D^ is given by 

1 



with a lattice spacing a, and U-^{x) = uji{x — jl). Here, /-matrix is defined to be hermitian, i.e., 
yji = 7^1- Thus, E) becomes an antihermitian operator, and Dirac eigenvalues are pure imaginary. 
We introduce the normalized Dirac eigenstate \n), which satisfies 

lD\n) = iX„\n), (2.2) 
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with A„ G R, and an eigenfunction i//,, (x) is expressed as 

Ynix) = {x\n), (2.3) 

which satisfies D\j/n = i^nWn- 

We introduce the operator formalism in lattice QCD [H, ^ 0], which is constructed from the 
link- variable operator . The link- variable operator is defined by the matrix element as 

{x\U^\y) = U^ix)5,+f,,y, (2.4) 

using the original link- variable U^{x). We define the Dirac-mode matrix element {n\Un \m) as 

{n\U^\m) =Y,{n\x){x\U^\x + fi){x + fi\m) = Y,¥l(.x)U^{x)Ym{x + fi), (2.5) 

X X 

with the Dirac eigenfunction Yn{x)- 

Considering the completeness relation ^„ \n){n\ = 1, any operator O can be expressed as 

d = ££|«)(«|d|m)(m|, (2.6) 

n m 

using the Dirac-mode basis. Note that this procedure is just the insertion of unity, and it is math- 
ematically correct. This expansion (2^) is the mathematical basis of the Dirac-mode expansion 
method [||, ^, ^. Next, we consider the Dirac-mode projection by introducing projection operator 
as 

for arbitrary set £/ of eigenmodes. For example, IR and UV mode-cut operators are given by 

PiR= £ \n){nl Ajv= £ \n){n\, (2.8) 
•^h|>Air |A„|<Auv 

with the IR/UV cutoff scale Air and Auv- 

Using the projection operator P, the Dirac-mode projected link-variable operator is given by 

Ui;=PU^P= £ £ \n){n\U^\m){m\. (2.9) 

We can investigate the Dirac-mode dependence of various kinds of quantities, e.g., the Wilson loop 
[^, ^, using the projected link- variable Uj^ instead of the original link- variable operator {/^. 

2.2 Polyakov loop operator and Dirac-mode projection 

Next, we formulate the Dirac-mode projected Polyakov loop. Here, we consider the periodic 
SU(3) lattice with the space-time volume V = x Nf and the lattice spacing a. In the operator 
formalism of lattice QCD, the Polyakov loop operator is given by 

1 ^' 1 

^^-^n^4=^t/f (2.10) 
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with the temporal link-variable operator U4. By the functional trace "Tr", the Polyakov loop oper- 
ator coincides with the standard definition as 

Trip = i-Tr{n(/4}= 3^ti-£(x,^in{/4lx,0 

i= 1 x.t i= 1 

■:^tTY,{x,t\U4\x,t + a) {x,t + a\U4\x,t + 2a) ■ ■ ■ {x,t + {Nt - \)a\U4\x,t) 



3V t 
1 

3V 



x.t 

^ trY^U4{x,t)U4ix,t + a)---U4{x,t + {Nt-l)a) = {Lp), (2.11) 

x,t 

where "tr" denotes the trace over SU(3) color index. 

We define the Dirac-mode projected Polyakov loop (Lp°^ ) as 

1 ^' 1 
^proj. ^ {Y\U^} = —Tv{PU4PU4P---PU4P} 

= 3^'^'' L {ni\U4\n2) {n2\U4\n3) ■ ■ ■ {nN,\U4\ni) . (2.12) 

In particular, the IR and the UV Dirac-mode projected Polyakov loop are denoted as 

{Lp)iK = ^tr £ {nr\U4\n2)---{nN,\U4\m), (2.13) 

|A„,.|>AiR 

(Lp)uv = ^tr £ {ni\U4\n2)---{nN,\U4\ni), (2.14) 

•^n,- |<AuV 

with the IR/UV eigenvalue cutoff Air and Auy. 

3. Lattice QCD calculation for Dirac-mode projected Polyakov loop 

In this section, we calculate the Dirac-mode projected Polyakov loop using SU(3) lattice QCD 



at the quenched level. We evaluate the full Dirac eigenmodes using LAPACK [12]. For actual 
calculation, we use the eigenmode basis of the Kogut-Susskind (KS) operator of 

= T- t W [Unmx+fi,y - 1/-^ (x)5,-_^,,] , (3.1) 

with 771 (x) = 1 and T]^(x) = (— l)^i+ '+^f -i > 2) in order to reduce the computational cost. The 



use of the KS-Dirac operator gives the same result as the original Dirac operator in Eq. ( p.l[ ) for 
the Polyakov loop. 

3.1 The confined phase 

First, we analyze the Polyakov loop in the confined phase below T^. Here, we use 6^ lattice 
with j8 = 5.6, which corresponds to the lattice spacing a ~ 0.25 fm and T = 1/ {Nto) ~ 0. 13 GeV 
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Figure 1: The spectral density p{X) of the Dirac operator on 6^ lattice with j8 = 5.6, i.e., a ~ 0.25 fm: (a) 
original spectral density, (b) IR-cut Pir(A) at Ajr = 0.5fl^', (c) UV-cut Puvi^) at Auv = 2.Qa^^. 
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Figure 2: The scatter plot of the Polyakov loop in the confined phase on 6^ lattice with j3 = 5.6, i.e., 
a — 0.25 fm and T = 1/ (Nta) ~ 0.2 GeV. (a) The original (no Dirac-mode cut) Polyakov loop, (b) The 
low-lying Dirac-mode cut at Air = 0.5a^'. (c) The higher Dirac-mode cut at Ayy = 2.0fl^'. 



Figure H shows the lattice QCD results for the Dirac spectral density p(A), and IR/UV-cut 
spectral density Pir(A) = p{X)G{\X\ - Air), Pi]y{X) = p(A)0(Auv — \M) with Air = 0.5a^^ and 
Auv = 2.0a^^. The total number of the KS Dirac-mode is x A/, x 3 = 3888, and both mode cuts 
correspond to removing about 400 modes. 

Figure ^ shows the scatter plot of the original Polyakov loop, low-lying Dirac-modes cut 
at Air = 0.5a^\ and the higher Dirac-modes cut at Auv = 2.0a^^, respectively. As shown in 
Fig. §^a), the Polyakov loop is almost zero, i.e., (Lp) ~ 0, which indicates the confined phase. 

Then, we consider low-lying Dirac-modes projection, which leads to the effective restoration 
of chiral symmetry breaking ^ [I^ In the presence of the IR cut Air, the quark condensate 
is given by 

(^^>IR = -^ I (3.2) 

At the IR cut parameter Air = 0.5a~^ ~ 0.4 GeV, only 2% of the quark condensate remains as 
{qq)m/ {qq) ^ 0.02 around the physical region m ~ 5 MeV |^]. However, as shown in Fig. §(b), 
the Polyakov loop (Lp)iR remains almost zero and Zs-center symmetry is unbroken, and these facts 
indicate that the system still remains in the confined phase, even without chiral symmetry breaking. 

In addition to the low-lying mode cut, we show the higher Dirac-modes cut in Fig. ^(c). In 
this case, the chiral condensation is almost unchanged, and the Polyakov loop remains almost zero. 
Therefore, the Polyakov loop is insensitive to both low-lying and higher Dirac eigenmodes. 
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Figure 3: The scatter plot of the Polyakov loop in the deconfined phase on 6^ x 4 lattice with j3 = 6.0, 
i.e., a = 0.10 fm and T = 1/ (Nta) ~ 0.5 GeV. (a) The original (no Dirac-mode cut) Polyakov loop, (b) The 
low-lying Dirac-mode cut at Air = 0.5fl^'. (c) The higher Dirac-mode cut at Auv = 2.0fl^'. 



3.2 The deconfined phase at high temperature 

Next, we investigate the Polyakov loop in the deconfined phase at high temperature. Here, we 
use 6^ X 4 lattice with j3 = 6.0, which corresponds to a = 0.10 fm and T = I /{Nta) ~ 0.5 GeV. 
The total number of the KS Dirac-mode is x M x 3 = 2592. 

We show the original Polyakov loop, typical low-lying mode cut at Air = 0.5a^^ and higher 
mode cut at Auv = 2.0a^^ in Fig. ^. These mode cuts correspond to removing about 200 eigen- 
modes. As shown in Fig. ^a), the Polyakov loop has non-zero expectation value (Lp) ^ 0, and 
shows the center group Z3 structure on the complex plane. These behaviors indicate the deconfined 
phase. 

As shown in Figs, ^(b) and (c), the Polyakov loop shows the characteristic behaviors in the 
deconfined phase even after removing low-lying or higher Dirac-modes. To be strict, the UV-cut 
Polyakov loop has a smaller absolute value than the IR-cut Polyakov loop, although the number 
of UV-cut modes is comparable to that of the IR-cut case. This suggests that contributions of the 
higher Dirac-modes are much larger than low-lying modes [Q] . However, apart from the normal- 
ization, both IR and UV cut Polyakov loops show the characteristic Z3 -pattern in the deconfined 
phase, and hence these Dirac-modes seem to be insensitive to the Polyakov loop properties. 

Figure ^ shows the Dirac spectral densities in confined and deconfined phases on 6^ x 4 lattice 
with j3 = 5.6 and 6.0, respectively. We also compare their low-lying spectral densities in Fig. ^c). 
In the deconfinement phase, the low-lying Dirac-modes are suppressed, which leads to the chiral 



restoration. The chiral condensate is also reduced by the IR cut of the Dirac-modes as in Eq. (|3.2|). 
On the other hand, there seems no clear correspondence between the Dirac spectral densities and 
the Polyakov loop. 

3.3 j3 -dependence of the Dirac-mode projected Polyakov loop 

Finally, we investigate j8 -dependence of the Dirac-mode projected Polyakov loop. Here, we 
adopt 6^ x 4 lattice with j3 = 5.4 ~ 6.0. 

Figure ^ is the absolute values of the Polyakov loop with typical Dirac-mode projections, 
and the original Polyakov loop data are also added for comparison. In this lattice volume, the 
deconfinement phase transition occurs around j3 = 5.6 ~ 5.7. As shown in Fig. S, both low-lying 
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Figure 4: The Dirac spectral densities in confined and deconfined phases, respectively (a) 6^ x 4 lattice with 
/3 = 5.6 in the confined phase, (b) 6^ x 4 lattice with j3 = 6.0 in the deconfined phase, (c) The comparison 
between confined and deconfined phases on low-lying spectral densities. 
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Figure 5: j3 -dependence of the absolute value of the Polyakov loop on 6^ x 4 lattice, (a) The IR Dirac-mode 
cut with Air = O.Sa^' and l.Oa^'. (b) The UV Dirac-mode cut with Auv = 2.0a"' and 1.7a"'. 



and higher Dirac-mode projected Polyakov loops show the similar j8 -dependence as the original 
data, apart from a normalization factor. This Dirac-mode insensitivity of the Polyakov loop is 
consistent with the results in the previous subsections. 

4. Summary 

In this paper, we have analyzed the Polyakov loop in terms of the Dirac eigenmodes using 
SU(3) lattice QCD. We have carefully removed relevant degrees of freedom for chiral symmetry 
breaking from the Polyakov loop. 

In the confined phase below T^, the Polyakov loop is almost zero, i.e., (Lp) ~ 0. By removing 
low-lying Dirac-modes, the chiral condensate {qq) is largely reduced. However, we have found 
that the Polyakov loop remains almost zero as (Lp)iR ~ even without low-lying Dirac-modes, 
and this fact indicates that the system still remains in the confined phase. We have also investigated 
contributions from higher Dirac-modes to the Polyakov loop, and have found no change of the 
Polyakov loop without higher Dirac-modes. Therefore, there seems no specific region of the Dirac 
eigenmodes essential for the Polyakov loop. 

We have also investigated the Polyakov loop in the deconfined phase at high temperature. In 
the deconfined phase, the Polyakov loop has non-zero expectation value, i.e., (Lp) ^ 0, which dis- 
tributes around Z3 elements in the complex plane. These characteristic behaviors also remain with- 
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out low-lying and higher Dirac eigenmodes. Therefore, the Polyakov loop and the Z3-symmetry 
behavior do not depend on low-lying and higher Dirac eigenmodes in both confinement and decon- 
finement phases. 

Here, we comment on the related studies about the correspondence between the Dirac eigen- 
modes and confinement. In the previous studies [|5[ ^, we investigated Dirac-mode dependence 
of the Wilson loop, and found that the confining potential survives without low-lying Dirac eigen- 
modes. The Graz group also reported that hadrons still remain as bound states without chiral sym- 



metry breaking by removing low-lying Dirac-modes [10, 1 1 1, which seems to suggest the existence 
of the confining force. 

These lattice QCD studies suggest that there is no direct relation between chiral symmetry 
breaking and confinement through the Dirac eigenmodes. For further investigation of correspon- 
dence between these phenomena, it is also interesting to analyze chiral symmetry breaking from 
the relevant eigenmodes of confinement [|l3[]. 
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